Excitation spectrum of doped two-leg ladders: A field theory analysis 
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We apply quantum field theory to study the excitation spectrum of doped two-leg ladders. It 
follows from our analysis that throughout most of the phase diagram the spectrum consists of 
degenerate quartets of kinks and anti-kinks and a multiplet of vector particles split according to the 
symmetry of the problem as 3 + 2 +1. This basic picture experiences corrections when one moves 
through the phase diagram. In some regions the splitting may become very small and in others 
it is so large that some multiplets are pushed in the continuum and become unstable. At second 
order transition lines masses of certain particles vanish. Very close to the first order transition 
line additional generations of particles emerge. Strong interactions in some sectors may generate 
additional bound states (like breathers) in the asymmetric charge sector. We briefly describe the 
properties of various correlation functions in different phases. 



PACS numbers: 71.10.Pm, ll.lOKk, 72.80.Sk 



I. INTRODUCTION 

The problem of ladder-like materials has attracted a 
lot of attention since the original paper by Dagotto and 
Ricei. Important experimental realizations of two-leg 
ladder systems include, for example, the famous 'tele- 
phone number' compound Sri4_. E Ca 3 ;Cu2404i which has 
a spin gap and exhibits a transition from Charge Density 
Wave (CDW) to superconducting (SC) state under the 
increase of Ca concentration2iiMiiii£iL2ai. Besides the ex- 
perimental relevance the study of this problem gives rise 
to many questions of rather general character. Even if 
ladder systems represent only the first step from a purely 
one-dimensional world into higher dimensions, this step 
introduces a lot of new interesting physics. In addition, 
the problem of doped spin ladders is just a particular 
case of a more general problem of multi-orbital quasi- 
one-dimensional models^, and an increase in the num- 
ber of orbitals dramatically increases the complexity of 
the lattice Hamiltonian, rising legitimate questions about 
universality 

In this paper we focus on the problem of doped two- 
leg ladders. Assuming that the spectral gaps are much 
smaller than the bandwidth (the applicability of these 
assumption to real systems is discussed in the end of the 
paper), we study the low-energy physics using the field 
theory approach. The intensive theoretical research con- 
ducted on this subject in the last eight years has estab- 
lished the following facts. 

A. Strong-weak tunneling duality. The form of the effec- 
tive action describing the low energy behavior of the sys- 
tem is independent (up to a simple operator transforma- 
tion) on whether one takes into account the inter-chain 
tunneling or considers just the inter-chain interactions. 
One arrives to this conclusion comparing effective actions 
deriv ed in the lim it of strong tunneling (see, for example, 
Ref.s ll lll2ll3ll4f) with the theories derived in the limit 
of weak tunnelingiSii&. The strong-weak tunneling dual- 
ity may be a boon for numerical calculations allowing to 



extract information about the excitation spectrum in one 
parameter range by performing actual calculations in the 
other. 

B. Superconductivity- CDW duality. The phase diagram 
includes the Tomonaga-Luttinger (TL) type phase as well 
as phases with spectral gaps. The only stable TL phase 
is the one where all modes are gapless. As soon as an 
attractive interaction appears in one channel it induces 
attraction in all other ones generating spectral gaps for 
all modes except the symmetric charge mode. The strong 
coupling phases are characterized by power law correla- 
tions for particular operators (order parameters). These 
phases are classified either as SC or CDW. For the two 
chain model there are two SC phases (s and d) and two 
CDW phases (also s and d, the latter phase also being 
known as Orbital Antiferromagnetii or Staggered FlvrsiSi 
phase) . Such classification is valid only for a weak repul- 
sion when the Luttinger parameter K c characterizing the 
gapless charge mode, is close to one. At K c < 1/2 the 
Wigner crystal phase with Akp density correlations com- 
petes with the SC phasesi2i2S. The low-energy Hamilto- 
nians in different sectors differ by the sign of certain cou- 
pling constants and transform to each other by canonical 
transformations of the fields. These transformations re- 
alize authomorphisms of the 0(6) groupie. 

C. Emergent attractive interactions. As follows both 
from the numerical^SiSi and analytic calculations based 
on the renormalization group (RG) analysis 1 2iiiiSi2i2^, 
the system at low energies may enter into a strong cou- 
pling regime even if all bare couplings are repulsive. In 
that case the system passes through an intermediate weak 
coupling regime. If E* is the energy at which the weak 
coupling is reached, than it can be shown that the strong 
coupling regime is achieved at energies of the order of 
E* 2 /A, where A is the bare ultraviolet cut-off. It is prob- 
ably safe to say that the emergent attraction leads to 
small gaps. 

It was first suggested in Ref. that if the bare cou- 
plings are not very large, the low energy sector of a 
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two-leg ladder is described by a universal Hamiltonian 
with a symmetry larger than the symmetry of the lattice 
model. This corresponds to the 0(6)xU(l) symmetry 
in the doped case and the 0(8) symmetry at half filling. 
The original suggestion was based on the observation that 
the RG flows of different coupling constants converge in 
the strong coupling limit to the same asymptotics cor- 
responding to a higher symmetry. Convergence of RG 
flows has also been found to occur in other systems^ and 
is frequently associated with Dynamical Symmetry En- 
largement (DSE). A careful discussion of DSE for some 
specific models can be found in Ref. |2(| This problem 
is addressed also here matching the RG analysis, valid 
at weak coupling, with methods valid at strong coupling. 
The conclusion is that in models where the number of 
particle flavors is not large, there are conditions when one 
may indeed expect a small splitting of particle multiplets. 
Such conditions exist in systems with weak backscatter- 
ing, like carbon nanotubes. 

The goal of this paper is to use a quantum field theory 
approach to describe the phase diagram and excitation 
spectrum of a doped two-leg Hubbard-type model. We 
combine RG with methods suitable to study the strong 
coupling limit, like X/N expansion^!, exact solutions and 
perturbation theory around specific integrable points in 
the parameter space of the effective field theory2£. We 
also discuss the conditions of validity of the field theory 
description. As said above, this description is valid when 
the correlation length is much larger then the ultra-violet 
(UV) cut-off. This regime is known as the scaling limit. 

The paper is organized as follows. In the next Sec- 
tion we introduce the effective low-energy description 
of doped two-leg ladders, the so-called generalized 0(6) 
Gross-Neveu model, and briefly analyze the RG equa- 
tions. In Sec. lIIII we discuss the possible strong coupling 
phases of the model and in Sec. IIVI we describe its sym- 
metries and identify the phase boundaries. The spectrum 
of generalized Gross-Neveu models is discussed in Sec. El 
using the large N approximation. Some specific points 
of the phase diagram are described by integrable mod- 
els with the 0(6) GN model being the most prominent. 
These integrable points are identified in Sec. IVII In the 
vicinity of these points we study the spectrum using a 
specific strong-coupling perturbation theory described in 
Sec. IVIII The spectrum close to the phase boundaries 
is summarized in Sec. I Villi Sec. IIXI is devoted to the 
study of a simplified 0(3)xO(3)-symmetric model where 
the RG equations can be integrated explicitly. In Sec. |3 
we briefly consider the structure of some correlation func- 
tions. We summarize the main results and discuss the ex- 
perimental relevance in the last section. The paper has 
several Appendixes. 



II. THE MODEL 

The model in its original formulation includes electron 
creation and annihilation operators c\ ; a , c^z^ labeled by 



chain indices, / = 1,2, and spin indices, a = +1 for spin 
up (t) and a = — 1 for spin down (|). For instance, the 
extended Hubbard model has the standard form 

H = -t|| (4,i, tr Ci + 1 < l << T + h - c ) + U ^2 n iM n U,l 

1,1, (J i,l 
-t X (c\ l l7 C h 2,a + h.c)j 

+V\\ ^ n i,l n i+i,l + V± E n M n *.2, (1) 
z,Z i,l 

where, as usual, the parameters fn, t± are the hopping 
matrix elements along and between the chains, U is the 
on-site repulsion and V\\ , V± the next-nearest neighbor 

interactions, with = c\ lrT Ci^ a and m t i = n ljiiT + 

riii^. The lattice Hamiltonian has a U(l) x SU(2) x 
Z2 symmetry. If necessary, one can include also exchange 
interactions. 

A. Low energy field theory 

As we have mentioned in the introduction, the effec- 
tive field theory, describing the low energy behavior of 
the Hamiltonian (JJJ, is largely independent on whether 
one considers strong interactions and weak inter-chain 
tunneling, as it was done, for example, in Ref.s Il5lll6l 
or diagonalizes the inter-chain hopping first and treats 
the bare interactions as weak, as it has been done in 
the majority of other papers. To be precise, there are 
two differences: (i) in one case the fields are labeled by 
chain indices and in the other by transverse wave vectors 
q = 0, 7r, (ii) the Hamiltonians of both sectors are related 
to each other by the particle-hole transformation. 

As far as the analysis of the phase diagram and the 
spectrum are concerned, it is more advantageous to 
have a low energy description in the Majorana fermion 
representation. The Majorana formulation is derived 
by bosonization and subsequent refermionization of the 
original Hamiltonian (the procedure was introduced in 
Ref. I2H see also lsnj) and, for th e spec ific model JD, can 
be found, for instance, in Ref.s Il4lll9l Away from half- 
filling, the resulting low energy field theory consists of 
two parts: one contains a decoupled symmetric charge 
mode $c ' and the other contains all other fields. The 
latter can be written in terms of 6 Majorana fermions as 

H = H + V, (2a) 

where Hq is the free part 
. 3 

a=l 
. 3 

- 55>.(&0«&-£0x&) (2b) 
0=1 
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and V describes the interaction between the right and 
the left moving Majorana fermions x R l anc ^ £fl l 

v = -g.+ixRxif - g P - [(^i) + (e R e L )] 2 
-2(x a Bxi) {9«-(&&) + [(e R e L ) + } 

-2ffc,«(£&i) + > (2c) 

(the summation over repeated indices is assumed). The 
symmetry of the continuum Hamiltonian is£7(l)x£/(l)x 
SU(2) x Z 2 , being somewhat higher than the symme- 
try of the lattice model. As we have already men- 
tioned, one U(\) field (the total charge mode $ c ) is 
decoupled and is not shown in the above Hamiltonian, 
while the second f (1) symmetry emerges only asymp- 
totically at lower energies. The mode $c is respon- 
sible for high conductivity along the ladders observed 
in Sri4^ x Ca a; Cu2404i 8 . The SU(2) triplet x° and the 
fermion £ 3 are made of bosonic fields of the spin sector. 
They reflect the SU(2) x Z 2 symmetry of the spin sec- 
tor and also appear in the low energy description of the 
two- leg Heisenberg ladder 29 . The situation here differs 
from the Heisenberg ladder since no explicit mass term 
is present in the Hamiltonian and the masses are gener- 
ated dynamically. The fermionic doublet £1.2 describes 

the asymmetric charge mode, which we will denote 0c \ 
The original observables are nonlocal in terms of the 
Majorana fermions (vector particles). The latter parti- 
cles, if remain stable, represent collective excitations of 
the system. We emphasize that the vector multiplet in 
is naturally split into sub-multiplets as 3 + 1 +2. 
The model © can also be represented as six critical Ising 
models coupled together by products of the energy den- 
sity operators. This representation is convenient because 
the original fermionic bilinears are local in terms of the 
Ising model order and disorder parameter fields. To clar- 
ify the symmetries and the structure of the model, we 
can rewrite it as 

H = H 0{3) [ X ;9a,+} + H 0(2) [e[r ) ;g p ^]+ H Ismg [f] 

+2i( XR xl){i9*-(&&) + 9c,st cos[/36 c -)]} 
+i2.g c , ss (^3 )cos[/30 (-) ] _ (3) 

Here Hi S i ng is the Hamiltonian of the critical Ising 
(CI) model and H Q (n) represents the 0{N) Gross- Neveu 
(GN) model^i that can be described in terms of N Ma- 
jorana fermions ifi^ L (a = 1, . . . , N) as 

h 0{ n)W,9n] = -^v N (nd x ^ R -rLd x r L )-9N(r R r L ) 2 - 

The spectrum of Hq(n) is massive for < 0, provided 
N > 2, while H (2) is always massless and is equivalent 
to the Gaussian model. 

Though the form of the Hamiltonians is fixed 

by the symmetry, estimates of the coupling constants are 
available only for weak interactions. The backscatter- 
ing interactions are weak in such systems as carbon nan- 
otubes which have the same symmetry as two-leg ladders, 



but not in two-leg ladders themselves. Therefore our phi- 
losophy will be the same as in the particle physics: we 
will express everything in terms of low energy parameters 
(mass gaps), and assume these gaps to be much smaller 
than the ultraviolet cut-off (whatever this cut-off is). As 
an additional simplification we will ignore the difference 
between the velocities, setting v s = v a = 1. 

B. Preliminary RG analysis 

The first step in our analysis is to establish the condi- 
tions under which the system scales to strong coupling. 
To do i t we use the single loop RG equations as obtained 
in Ref . Illll2ll6l2lt 

9 P - = -3#c, s t - 9c, ss (5a) 

9c,ss = -9 P ~9c,ss - ^9a-9c,st (5b) 

9<r- = -2ga+ga- - 2g c ^ t g c ,ss (5c) 
9c,st = -{g P - + 2g a +)g c ,st - ga-gc,ss (5d) 
fU+ = -gl+ - gl- - 2g 2 c st . (5e) 

Here the dot corresponds to a derivative in RG time t = 
(47rw) _1 ln(A/e). Thus small energies in our notations 
correspond to large t. 

Further in the text we will study the RG equations 
close to the 0(6)-symmetric point and, in more detail, for 
a simplified 0(3)xO(3)-symmetric model (cf. Sec. 11X11 . 
In the latter case Eqs.© can be solved analytically. 

It is possible to show that there are two areas of sta- 
bility. 

• The weak coupling area, where the repulsion domi- 
nates. All couplings scale to zero. The resulting 
phase is the 0(6)xU(l) TL liquid perturbed by 
marginally irrelevant perturbations. We found the 
exact boundaries of this area for the 0(3) x 0(3) 
model, but qualitatively one can say that the TL 
liquid appears as a fixed point when the diagonal in- 
teractions are repulsive and exceed the off-diagonal 
terms: g a+ ~ g p ,- > \g„-\, \9c,ss\, \g c ,st\- 

• The strong coupling area. In the RG sense this area 
is the basin of attraction of the 0(6)-symmetric 
point -g a+ = -g p ._ = \g a _\ = \g c . ss \ = \g c . st \ > 0. 
In order to get to strong coupling it is enough 
to have just one attractive diagonal interaction. 
Moreover and very intriguingly, one may even start 
from all interactions being repulsive. The sys- 
tem will scale to strong coupling anyway provided 
the off-diagonal interactions are sufficiently strong 
(this would correspond to emergent attraction de- 
scribed in Introduction). This follows from the 
stability analysis of the unstable symmetric line 

ga+ = 9 P - = \9(t-\ = \g c ,ss\ = \g c ,st\ > 0. 

The RG equations determine an overall scale at which 
the strong coupling is achieved: 

M = A/[. 9a (0)] = A/(. 9l (0), d, ...C Q -i). (6) 
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Each set of initial conditions sets the system on a par- 
ticular RG trajectory defined by its invariants Ci(i — 
1, ...Q — 1), where Q is the number of the coupling con- 
stants. For a given set of Cj's the scaling limit exists if 
for an arbitrary large A one can choose a starting point 
on the trajectory such that M remains constant. It turns 
out that for trajectories with emergent attraction this is 
not possible. Therefore the space of C,'s where the scal- 
ing limit is defined is a subspace of (Q — l)-dimensional 
space (scaling subspace). Inside of this space the masses 
may depend on C,-'s and this dependence is not deter- 
mined perturbatively. In view of this it is clear that the 
existence of "strong coupling point" or DSE requires in- 
dependence (or at least weak dependence) of the mass 
spectrum on RG invariants. Otherwise strongly corre- 
lated phases are not points, but have their own geogra- 
phy with the excitation spectrum changing throughout 
the phase. We will continue this discussion in Sections 
IVII CI and IIXI where we consider mass variations related 
to small deviations from the 0(6) and 0(3) x 0(3) sym- 
metries respectively. 



III. STRONG COUPLING PHASES AND 
ORDER PARAMETERS 



To explore the possible phases of the model it is more 
convenient to use the mixed representation @. Here 
the Hamiltonian is expressed in terms of bosonic fields 
associated to the charge modes and Ising fields a a 
(fjL a ) with a — 0, . . . , 3 associated to the spin modes £ 3 
and x° respectively. The phases of the system are in 
one-to-one correspondence with the fields vacua. These 
vacua are determined by the signs of the renormalized 
couplings (g a -,g c . s t,g c .s S ) in the strong coupling regime 
(recall that these signs may have nothing to do with signs 
of the bare coupling constants). The analysis of the phase 
diagram has already been conducted and here we repeat 
many results obtained in Ref. Il3lll4l ll6 20. To keep con- 
tact with these works we use a strong tunneling approach 
introducing bonding (p = 1) and anti-bonding (p = —1) 
operators 



0(n)p, ( 



Cn,l,a + P C W ,2,<T 

V2 



(7) 



associated to transverse wave vectors q — 0, ir respec- 
tively. The bosonization notations are given in Appendix 
A and the relationship with the weak tunneling approach 
is discussed at the end of the section. 

When the forward scattering in the symmetric charge 
channel is not strong, there are four possible phases: su- 
perconducting s- (SCs) and d-wave (SCd), charge den- 
sity wave (CDW), and, what is now frequently called, 
d-wave CDW (CDWd). Its order parameter is of Orbital 
Antiferromagnet 17 -Staggered Flu»i& type. In terms of 
the lattice operators the order parameters of the four 



phases have the form 

AsCs(«) = ^ C ( n )p,<r C { n )p,-<r 
p=±l 

P =± 

&scd{n) = ^2 sm(Trp/2)C(n) p , a C(n) Pt . 



p=±i 



y"l P [Rp^L p .i + Lp^R 



p,1\ 



(8a) 



(8b) 



P =± 



A CDW (n)= J2 CHn) Pta C(n)^ a e 2k r ina ° 

p=±l,<7 

= Yl K* L -p,° ( 8c ) 

p— ±1.(7 

A C DWd(n) = sin(7rp/2)O t (n) PiCT O(7i)_ P:ff e 2 ^ i " Q0 

p=±l,a 

= J2 p R t° L - P ^- ( 8d ) 

p=±l,cr 

Here we introduced the right and the left moving com- 
ponents of the lattice fermion operators: 

C p , a {n) = e- lk F )x R p ^(x) + e lk ? x L p ^(x), (x = na )(9) 
(p) 

where kp are Fermi vectors of the different bands and 

k F = (fc^ +) + k i f ) )/2 = tt(1 - d)/2a , with S being the 
doping and ao is the lattice constant. 

The SCd phase is found in the region (+, +, +) and is 
characterized by 

e£-> =0> o ) ^0(a = 0,l,2,3) or 
6<-> = V5F/2, </J a ) ^0(a = 0,1, 2, 3) (10) 

where (...) = (0| . . . |0) is the vacuum expectation value 
(VEV). This phase has power law correlations (quasi- 
long-range order) in the d-wave Cooper channel 18b(l . 
The order parameter in the continuum limit is 

A scd ~e^ +) (11) 
x |cos[\/7rO^ :) ]o-icr2(T30'o - isin[v^O£~^]jUiM2M3Mo} • 

For (-,+,-) one finds the CDWd phase, characterized 
by 

e£-> = 0, (a a ) i= (a = l,2,3),/i Q ^0 or 
O(-) = v^/2, (/ia) ^0 (a = 1, 2, 3),a ^ 0. (12) 

One gets to this phase from the first one through a Z 2 
QCP. In the continuum (|8d|) becomes 

AcDWd ~ e-^ +) (13) 
x |isin[V7r©i _) ]MiA*2M30'o - cos^S^T Vl^CTaMo} • 



From the CDWd going through a U(l) QCP one 
reaches the CDW phase characterized by (-,-,+) and 

0(~) = (a a ) + 0(a = 1, 2, 3), M o + or 

6(") = 0, ( Mo ) ± 0(a = 1, 2, 3), (7 ^ 0. (14) 

The order parameter in the continuum limit takes the 
form 

Acdw ~ e^ +) (15) 
x |cos[-\/7r6i~^]/iiM2A t 30-o - sin[y/Tr@^]aia 2 (T 3 fi \ • 

Finally, the SCs phase is dominant in the region (+,- 
,-) and in characterized by 

0(-) = V5F/2, (<r ) ^ 0(a = 0, 1, 2, 3) or 

e£-> = 0, < Ma ) ^ (a = 0,1, 2, 3). (16) 

Its order parameter in the continuum limit is 

Asc.^^ (17) 
x |-isin[\/7re^ :i ](TiCT2cr30'o + cos[Vti : B^Vi^2M3 ) "o} • 

This last phase can be reached both from the SCd phase 
through a U(l) QCP or from the CDW phase through 
a Z2 critical point. All possible phase transitions are 
schematically shown in Fig. ^ and can be summarized 
as follows: Z 2 QCP between SCd and CDWd, SCs and 
CDW; U(l) QCP between SCd and SCs, CDW and 
CDWd and 1st order transition between SCs and CDWd, 
and between SCd and CDW. More details about the 
phase boundaries are provided in the next Section. 

The scaling dimension of both CDWs phases is 
d-CDW — K c /4, while the one of SC phases is dsc = 
1/AK C . The obtained bosonized expressions for the 
CDW order parameters differ from the ones obtained in 
Ref.slHi^ In our treatment of the order parameters we 
used Ising model notations. This is done because Ising 
model order and disorder parameters a and \x go natu- 
rally with Majorana fermions. 

The density operator also contains a Akp — 2{kp + 
kp ' ) oscillatory piece 

p{Ak F ;q) = iexp[i^$(+)] [A q ( x a RX a L ) + B 9 &&] ,(18) 

where q — 0, ir is the transverse momentum. Notice that 
the quasi-long- range order occurs only at the wave vec- 
tor corresponding to the overall particle density. Cor- 
relations with 4fc^ l ~' , 4k F ^ turn out to be exponentially 
suppressed. The amplitudes Aq , Bq vanish for the non- 
interacting system and for weak interactions they are of 
the order of U/cf- The impo rtance of the 4k p correla- 
tions was pointed out in Ref.s ll9l20i Please notice that 
since the Majorana bilinears always have nonzero expec- 
tation values as soon as the spectrum is gapful, the am- 
plitude of the 4fcp wave is finite except in the Tomonaga- 
Luttinger phase and therefore critical Akp density fluc- 
tuations are always present. As a consequence the 4/c^ 
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Isin g SCd U(l) 




FIG. 1: The two-dimensional projection of the phase diagram 
for K c > 1/2. The thick lines represent first order phase 
transitions. At K c < 1/2 the SC phases are replaced by the 
Wigner crystal. 



order parameter, having the scaling dimension d± = K c , 
competes with the SC ones. At K c < 1/2 the scaling 
dimension of l|18l) becomes smaller than the scaling di- 
mension of the SC order parameters and SC phases on 
the phase diagram are replaced by the Akp CDW which 
we call Wigner crystal^. It is possible that this mecha- 
nism is responsible for charge ordering in the telephone 
number compound observed in Ref. 7 6 9. As far as the 
CDW phases are concerned, the 2kp oscillations are al- 
ways more relevant. However, an establishment of a true 
long range 2kp CDW order in a three-dimensional array 
of ladders will also lead to condensation of the 4/cp os- 
cillations. At special value of 4k p — 2tt (the 1/4 filled 
band) operator (|18(l couples to the lattice. At this point 

the charge field couples to the other fields and the 
corresponding sector acquires a gap. As follows from nu- 
merical calculations done in Ref. l33l in two-leg ladders 
this may occur in a broad range of parameters. 

As far as higher harmonics of the electron density with 
the wave vectors 2nfcj? are concerned, their scaling di- 
mensions grow as n 2 K c /4 and they quickly become irrel- 
evant. For this reason X-ray scattering from such Wigner 
crystal shows sinusoidal oscillation of the electron den- 
sity. This picture holds when the dominant interactions 
are smaller than the one-dimensional bandwidth so that 
the continuous description can be used. 

The order parameters in the model with weak inter- 
chain tunneling (smaller than the gaps) can be obtained 
from the above order parameters by the chiral particle- 
hole transformation: 



Ri,„ -> <rB£_ a (19) 

Recall that in the transformed theory I becomes a chain 
number and $ and O fields in the charge sector inter- 
change. 
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IV. DISCRETE SYMMETRIES AND PHASE 
BOUNDARIES 

Apart from the continuous symmetries mentioned 
above, the Hamiltonian (J2J possesses a set of discrete 
symmetries. These symmetries establish a one-to-one 
correspondence between excitation spectra of the differ- 
ent phases and represent authomorphisms of the 0(6) 
groupie. To see this we rewrite the interaction (|2*c)) in 
terms of S0i(6) Kac-Moody currents J a — ipRT a ipR and 
J a = ipLT°"ipL, where ipR^L represent all six Majorana 
fermions and r° are generators of the 0(6) symmetry 
group 

V = GabJaJb. (20) 

The discrete transformations which leave the Hamilto- 
nian invariant correspond to a change of sign of some 
chiral currents and the corresponding coupling constants, 
leaving the currents of opposite chirality unchanged 

Ja J(i: Ja Ja ■ (21) 

Since the transformations (|21|l must preserve the commu- 
tation relations of the currents, they correspond to autho- 
morphisms of the corresponding group (the 0(6) one in 
the given case) . The above authomorphisms do not affect 
the spectrum, but change the order parameters (recall 
that the latter ones are nonlocal in the Majoranas) and 
therefore establish a duality between different phases. In 
terms of Majorana fermions the authomorphisms corre- 
sponds to a sign change of some Majorana fermions in 
one chiral sector. Therefore as far as the spectrum is 
concerned, one can get a complete picture by studying it 
in a given phase, for instance, in SCd. All other phases 
can be obtained using the duality transformations. In 
the given case we have three dualities 

• £l * £l ' 9a, — i9c,ss * ~9cr,— i~9c,ss 

• xl^-xh 

9c : st-> {JciSS > Qc,st-> 9c, ss 

• Xl -> 9*,-,9c,st -> -9a-,~9c,sf 

Let us now consider in more details the boundaries 
between the different phases. Since the spectra of differ- 
ent phases are the same, we can take SCd phase as an 
example. 

The boundary corresponding to a Z 2 QCP separating 
SCd phase from OAF is a surface in the coupling con- 
stant space on which the Majorana fermion £ 3 decouples 
from the rest and becomes massless. This surface is de- 
termined by the condition 

g„-{{x a R xi)) +ffc, M ([«ki) + (eiei)]) = (22) 

The criticality is not violated by fluctuations around i|22|) 
since they cannot generate any relevant operators. 



Analogously, if (xrXl) an d (£j?££) take expectation 
values and 

9c, s t{(XRXl))+9c, S s((ZUl))=0, (23) 

the Majoranas £\£ 2 decouple. This is the U(l) QCP 
between SCd and SCs. The fluctuations generate a per- 
turbation cos[\/167r/30c '] which may become relevant if 
f3 2 < 1/4. This, however, corresponds to a very strong 
9 P ,— - Then this critical line may become a first order 
transition. 
Finally, if 

9*,-{(&&)) + (£&£)]) - 0, (24) 

with different from zero, the Majorana triplet \ a 

is decoupled from £ fermions. At g a ^ + < both models 
are massive, which means that this is a first order transi- 
tion. The fluctuations shift g a> + further to the negative 
side. If the effective coupling is positive, however, this be- 
comes a SU {2)2 quantum critical point. The first order 
line separates SC and CDW phases with different point 
group symmetries (for instance, SCd and CDW). 

For the analysis of the spectrum it is convenient to 
generalize the model © to N species of fermions, and 
rewrite the interaction (|2*c|l as 

v = -\x a lab x b ■ x a = ir R r L - (25) 

The vicinity of the symmetric line of the phase diagram 
is described by a weakly anisotropic O(N) GN model 

lab = -Tr(gO + $9ab) , ^9a,bl ' 90 < 1- (26) 

In an analogous way, at the phase boundaries the N x N 
matrix 7 is split into two N\ x N\ and N2 x N2 blocks 
7 1 and 7 2 respectively 

7° = (t 7 ° 2 ) ; lib = 9i , lib = 92, (27) 

and the vicinity of the phase boundaries is characterized 
by 

la,b = la,b + S 9a,b , £ffa,&/01,2 < 1- (28) 

For the case at hand we have N = 6, and 

N-l = 5, N 2 = 1, for Z 2 (29) 
N t = 4 , N 2 = 2, for 17(1) (30) 
N t =3, N 2 = 3, for the first order (31) 

boundaries respectively. 

V. SEMICLASSICAL ANALYSIS 

In this Section we provide some general arguments con- 
cerning the spectra of generalized O(N) GN models l|25t 
I28fl based on the semiclassical analysis valid in the limit of 
large NSL. A more intuitive analysis valid in the same ap- 
proximation but for even N, M is discussed in Appendix 

m 
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A. Effect of small anisotropy on the vector 
particles 

Let us consider first the generalization (|26|l of the GN 
model with a total number of fermion species N large. 
If the quadratic form (|25|l is negative defined, then it is 
possible to perform the Hubbard-Stratonovich transfor- 
mation introducing an auxiliary field A a 



V = ^A a ( 1 - 1 ) ab A b -A a X a , 



(32) 



where A is related to X by the equation of motion 

A a = lab X b . (33) 

Integrating out the fermions one finds the effective po- 
tential for A° 

1 N / A 2 \ 

V eff [A] = -A^k^-H-A 2 (log—- lj , (34) 

that has minima at A a = ±m a given by the saddle point 
equation 



a = 77- y^7 a &m b ln(A/|m b |). 

Z7T c — ' 



(35) 



From l|33(l and (|34|l it follows that the fermion bilinears 
acquire vacuum expectation values. 

Besides the trivial zero-energy solution, the classical 
static equation of motion for the field A in the potential 
(|34|l also has finite-energy, kink-like solutions, A c (x), in- 
terpolating between different minima 1)35(1 . When quan- 
tizing the theorji^iiS one finds that it possesses fermionic 
excitations with masses m , associated to the ordinary 
vacuum sector, and kink excitations, associated with con- 
figurations A c (x). In fact, as it was shown by Jakiw and 
Rebb: 36 , fermions interacting with topological kinks like 
in ll-'il'l possess a single normalizable zero energy mode, 
-00, in addition to the finite energy solutions ip n . The cru- 
cial point is that, while finite energy solutions are com- 
plex, the zero energy one is real and non-degenerate. The 
semiclassical expansion is then given by 3 ^^ 



-itE n 



- h.c. 



(36) 



where ip° is the zero energy solution of the Dirac equation 
and, assuming that A c (+oo) > 0, it has the form 



^ = rr 1 / 2 



exp 



(37) 



with f2 -1 ' 2 being the normalization factor. Since tp a are 
Majorana fermions, the operators "f a compose a Clifford 
algebra {70,7b} = <5 ao . Therefore these matrices realize 
spinor representations of the O(N) group. It is important 
to stress that the zero mode solutions exist for any con- 
figuration A(x) which has asymptotics of different sign (a 



kink). The effect of small anisotropy on the kink masses 
is difficult to study in the semiclassical limit; this prob- 
lem will be discussed in Sec. IVII CI using other methods. 
We now consider in detail the effect on fermions. 

The fermion masses are given by the saddle point 
equations l|35|) . The solutions of these equations de- 
pend on the bare couplings; for Sg a b = one recovers 
the usual form for the large N limit GN mass, namely 
m a = Mo = Ae~ 2 */ 90 . If 5g ab is non vanishing but 
\Sg a b\/go << 1, the solution is 



8m a 
Mo 



= % 



c b^c 



(38) 



From this analysis it is clear that the vector multiplet 
is split, the splitting being proportional to the splitting 
of the bare coupling constants and survives in the limit 
l/N -> provided TV -1 £ Sg ab ^ 0. 

Let us now imagine that some of the eigenvalues of 
7 are zero. The limiting case is when the matrix 7 is 
proportional to the projector: 



lab 



(39) 



Then Ea. l|35|l can be solved explicitly. The Hubbard- 
Stratonovich transformation yields 

£ = ^A 2 + (-^V - ie a Ar^ (40) 

and the masses of the vector particles are given by 

m a = e a (A) (41) 

where 

(A) = Aexp (-^ + A\ , A^N-'J2 e i WVM) 

The N — ► 00 limit is defined as follows: 

= N f dep(e), / dee 2 p(e) = 1, 

a J J 

A = J dep(e)e 2 ln(l/|e|) 
So N — 1 components of vector e are RG invariants. 

B. Kink confinement 

Let us now consider the opposite limit, when two 
0{Ni) and Q(N 2 ) GN m odels (N h2 > 2) are weakly 
coupled as in Eq.s I|25I28() . We discuss here only the evo- 
lution of the kink excitations since the analysis for the 
fermions is similar to the one carried out in the previous 
section. 
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In the absence of coupling, g a b = 0, each GN model 
has its own kinks, associated to two classical solutions, 
Ai^a:), and interpolating between the minima ±mi,„ 
and ±m2, a respectively. The effect of the interaction is 
to lift the degeneracy between the minima introducing a 
confining potential between the kinks. As a consequence 
the original kink solutions of the decoupled theories be- 
come unstable and disappear from the spectrum. 

To see this let us consider, for instance, the effect of the 
interaction on a kink-antikink configuration, where both 
particles belong to the same GN theory, while the other 
theory is at the minimum, for example A 2 (x) = m 2 . 
The kink-antikink configuration is such that Ai(x) takes 
value —mi from minus spatial infinity to a point x%(t) 
where it switches to +mi; it keeps this value until X2{t) 
where it switches back to —mi- In presence of the per- 
turbation, the kink-antikink state acquires an additional 
energy 



U(xi,x 2 ) = Sg a ,b\xi - x 2 \mim 2 . 



(42) 



This results in a confining potential for the two 0(N\) 
kinks (this universal form of the confining potential is 
valid only when the confining radius is much greater than 
the size of the lightest kink). For Ni > 4 this potential 
exists on top of the attractive potential already present 
in the 0(N\) GN model. Therefore the lowest bound 
states have the same symmetry as in the O(Ni) model, 
i.e. they transform according to the vector representa- 
tion of the group. Of course, if the confinement radius is 
much greater than the inverse kink's mass, the confine- 
ment potential contains many other bound states (we re- 
fer the reader to Sec. lVII Dl for a more detailed analysis). 
These states do not exist in the 0{N\ + N 2 ) GN model. 
Their masses are grouped around the mass of the 0(N\) 
vector particle m v i or, if N\ < 4, around 2M 1: where Mi 
is the kink's mass in the 0(N\) GN model. Repeating 
the same argument for a kink-antikink configuration be- 
longing to the 0(N 2 ) GM model one obtains the states 
grouped around the mass of the 0{N 2 ) vector particle 
m v2 . 

The same confining phenomenon happens also between 
two kinks that belong to different GN models, in this case 
the lowest multiplet of the confined states transforms ac- 
cording to the spinor representation of 0(Ni +N 2 ). This 
follows from the fact that such states correspond to a si- 
multaneous change of sign of Ai (x) and A2 (x) , therefore 
they are kinks of the 0(Ni + N 2 ) GN model. Since the 
only interaction between kinks of the different models is 
proportional to <5<? a ,b, the masses of these particles group 
around the sum of the kinks masses Mi+M 2 . The details 
of the structure of the excitations induced by the confin- 
ing potential is discussed in Sec. lVIII One should keep in 
mind that some multiple kink representations may sur- 
vive (see Appendix lB|l . 



C. Coupling massless and massive GN models 

The above analysis fails if the matrix 7 in Eq. (12511 has 
negative eigenvalues. If all eigenvalues are negative, the 
analysis of the RG equations shows that all interactions 
scale to zero leaving the theory massless. The situation 
when the form has no definite sign (at least, this is the 
criterion in the large N limit) corresponds to the area of 
emergent attraction. The large N analysis is still appli- 
cable here, but in a modified form. Let us consider the 
0(Ni) x 0(N 2 ) model with Ni/N 2 finite in a situation 
such that the interaction among the first N\ particles is 
attractive and the interaction among the other particles 
is repulsive. We take a simple form of the interaction 
to keep the analysis more clear. Then we can do the 
Hubbard-Stratonovich transformation in the attractive 
sector, so that after some algebraic manipulations we ob- 
tain (g+,g- > 0): 



- — A + iA £ a £ a xbXb 

2.9+ V 9+ 

AV 1 (J +.9-) (iXbXbf. 



9-c- , 2 
Wi ^XbXb) 



(43) 



Replacing A with a constant Mq m m^/g + where is 
a mass of £ fermions, we obtain the following Lagrangian 
for x fermions: 



£ = -jXal^Q^Xa + iMoXaXa ~ 9eff(XaXa) 



(44) 



where g e ff = N^ 1 {g\/ g + + .9-). Summing the diagrams 
with a maximal number of loops, we obtain the following 
expression for the mass ratio: 



mc 



9x_ 

9+ 



l + ^r(9 2 x/9 2 + + 9-/9+) 



(45) 



Thus also in this case the vector multiplet remains split. 

Let us now consider kinks. When one of the GN mod- 
els is massless the potential <|42ll can no longer be used. 
Like above, we can do the Hubbard-Stratonovich trans- 
formation in the massive sector (cf. Eq. 143(1 , and study 
the problem of the two types of fermions in a background 
bosonic field. For both fermionic modes the solution is 
again described by Eq. (|36I37|I . As was discussed in Sec- 
tion |^2J the entire kink is a bound state of the scalar 
field A (a;) and the fermions. The fermionic zero modes 
now realize spinor representations of the 0(N± + N 2 ) 
group and then in principle can be seen as bound states 
of massive kinks of the 0{N{] model and massless one of 
the 0(N 2 ). 

The summary of these semiclassical arguments is the 
following. Weak interaction between GN models gener- 
ate (i) degenerate kink multiplets transforming according 
to spinor representations of the 0(Ni + N 2 ) groups and 
(ii) vector particles with different masses transforming 
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according to the vector representations of the 0(N±) and 
0(N 2 ) groups. On top of it there may be many other 
multiplets with masses laying between (Mi + M 2 ), m„i, 
m V 2 and 2 (Mi + M2). These multiplets cross into the 
continuum and become progressively unstable once the 
coupling between two models increases. 



VI. INTEGRABLE POINTS 

At some specific points in the parameter space the 
Hamiltonian is intcgrable. One can identify the fol- 
lowing integrable models, that describe either the most 
symmetric points or the phase boundaries. 

0(6) GN model. As already noticed, if all coupling con- 
stants are equal Eq. J5J has an extended 0(6) symmetry 
and becomes the 0(6)-symmetric GN model Q. The 
0(N) GN model is integrable for any iV 34 i 38 i 39 . For even 
N the scattering theory is relatively simple^, while for N 
odd there are significant complications (see Ref . f4l|) . For 
g < and N > 2 the spectrum is massive and consists of 
kinks and antikinks, which transform according to the ir- 
reducible spinor representations of the O(N) group. For 
N even kinks and antikinks correspond to different rep- 
resentations, for N odd there is just one representation. 
The two-loop RG gives the following expression for the 
mass: 

M = A 3 1 /^ 2 )exph2 7 r/(A^- 2)g]. (46) 

For N > 4, there are also fermion particles and their 
bound states. These fermions correspond to the original 
Majorana fermions and transform according to the vector 
representation of the group. The masses of the fermions 
and their bound states are given by 

M a = 2Msm(jf^, a=l,... ) int(^). (47) 

For N even, the total number of kinks is 2 N / i ( reflect- 
ing the fact that there are two 2 JV / 2 ~ 1 -dimensional irre- 
ducible spinor representations of the algebra SO(N)SL). 
For N odd the total number of kinks is 2^ N+1 ^ 2 although 
there are major subtleties associated to multi-particle 
states^. The number of fermions is always N. In par- 
ticular, for N = 6 the spectrum consists of two spinor 
multiplets of mass M and one six-fold degenerate vec- 
tor multiplet with mass \[2M. No fermion bound states 
are present. An intuitive picture of the two types of ex- 
citations is provided in Appendix [B] for N even. It is 
also good to remember that the spinor representations of 
0(6) are isomorphic to fundamental representations of 
the SU(4) group; their quantum numbers correspond to 
spin a = ±1/2 and transverse momentum p — 0,7r (or 
chain index for theories with weak inter-chain tunneling) . 

0(3) x (U{1) x Z 2 ) model. For = g c<st = two 
groups of Majoranas decouple from each other; one is 
described by the 0(3) GN model, the other one by the 



anisotropic 0(3) GN (U(l) x Z2-symmetric) . Both mod- 
els are integrable fsee l45l46|) . and the excitation spectrum 
contains only kinks. 

(3+1) model. At the U(l) QCP, described earlier in 
Section III, two Majorana fermions £1,2 decouple and re- 
main massless. The remaining massive theory with the 
0(3) x Z2 symmetry is related to an integrable model 
solved by TsvelikiS, and Andrei, Jerea^ 3 .. (Some correla- 
tion functions for this model were calculated in Ref. IH .) 
The spectrum of this integrable model contains an SU (2) 
kink doublet with mass M and a light Majorana fermion 
of mass mo ~ M(\ — A'). There are no stable vector 
particles except of this singlet fermion who is, however, 
is not a bound state of kinks (see discussion at the end 
of the next paragraph). The S-matrix can be found in 
App.0 

(5+1) model. At the Ising QCP one Majorana decou- 
ples. For g a .+ = g P - = g c .st the massive model is 0(5) 
symmetric and integrable. Outside of QCP, provided the 
0(5) symmetry is maintained (g c , S s — 9a,-), it is possi- 
ble to construct a factorized scattering (integrable) the- 
ory with the same symmetry, 0(5) x Z2, and UV central 
charge. The exact solution is described in Appendix IU1 
the spectrum consists of two vector particles (quintet and 
singlet of Majorana fermions) with different masses mo 
and V3M and a quartet of kinks = antikinks with mass 
M realizing the spinor representation of the 0(5) group. 
This theory certainly deserves further analysis. There 
are two serious qualitative differences between the scat- 
tering theories for the (5+1) model and the 0(6) GN 
model (similar differences are found between the (3+1) 
model above and the 0(4)). First, as it follows from the 
group theory, in the former case the kinks and antikinks 
are the same particles and in the latter case they belong 
to different representations. Second, in the (5+1) theory 
the singlet fermion does not appear as a bound state of 
kinks as all vector particles do in the 0(6) GN model. 
As a consequence the singlet mass is not related to the 
kink's mass. One would imagine however, that with in- 
crease of the coupling between the Majorana singlet and 
the other ones the model approaches the 0(6) GN model. 
At present we do not have a complete picture of recon- 
ciliation of these two models. The most likely solution to 
this puzzle is that the two theories are equivalent only in 
the region where the singlet mass is small. 



VII. PERTURBING AROUND INTEGRABLE 
POINTS 

In this Section we discuss how to calculate the spec- 
trum using perturbation theory around integrable points. 
Since integrable points describe the phase boundaries and 
highly symmetric points, we will be able to obtain some 
additional information about the excitations of the model 
in these regions. 
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A. Form-factor perturbation theory 

Let us first recall some essential features of integrable 
models that we need in order to construct perturbation 
theory. Massive integrable models are characterized by a 
simplified on-shell dynamics which is encoded into a set 
of elastic and factorized scattering amplitudes of their 
massive particles. A convenient formalism for the de- 
scription of a dilute gas of particles with factorized scat- 
tering can be constructed in terms of the creation and 
annihilation operators, A\{6), A a (0), that satisfy the 
Zamolodchikov-Faddeev (ZF) algebra. (For an introduc- 
tion on these concepts see for instance Ref.s or the 
introductory chapters of Smirnov's book in 1471) Here the 
rapidity 9 parameterizes the relativistic dispersion rela- 
tion 

e a {9) = M a cosh 9, Pa {9) = M a sinh 9 (48) 

and a is an isotopic index. The ZF operators are the log- 
ical extension of the algebra of free fermions or bosons to 
the case of interacting particles with factorizable scatter- 
ing, where the interaction is completely characterized by 
the two-particle S-matrix, S at b(dij)**. As usual, multi- 
particle states are obtained acting with strings of creating 
operators on the vacuum 

\9 1 ,...,9 n ) ai _ an =A{ i (9 1 )...A{ n (9 n )\0). (49) 

The single particle states can be thought as generated 
by an operator, (p a (x), such that (0\(p a \9) a 7^ 0. For the 
O(N) GN model ZF operators include fermion, Af i (and 
their bound states) and kink, Ak t , operators. Though 
for the O(N) GN model the S-matrix is not diagonal, 
this will not be essential for what follows, so we prefer to 
describe the methods for diagonal S-matrices since the 
notations are less cumbersome. 

In what follows we will need a definition of non- 
locality. Recall that two operators, 0\ and O2, are 
said to be mutually non-local if the Euclidean correla- 
tor (. . . 0i(x)02(O) • • .) is not a single valued function 
of x. In particular, if we introduce complex variables 
z = x 1 + i x 2 and z = x 1 — i x 2 , and under analytic con- 
tinuation z — » ze 27rl , z — » ze _27rl , the correlator acquires 
only a phase, 2^71^, the two operators are said to be 
semi-local. The non-local operators that we will consider 
in the following are of this type. The great usefulness 
of this definition is that the index 7^2 can be calculated 
in the ultraviolet, where all correlation functions have a 
simple power law form. As we will see, the effect of a 
perturbation on the spectrum of a theory will crucially 
depend on the locality properties of the perturbing oper- 
ator. 

Let us imagine now that we perturb an integrable the- 
ory, Hi n t, with a non- integrable perturbation 

H = H int - g f dx V(x), (50) 

where ^(x) is a scalar field. The variation of the spec- 
trum of the theory can be studied perturbatively in g, 



in the same spirit as standard quantum mechanics (QM) 
perturbation theory, taking advantage of the fact that in 
integrable models matrix elements of perturbing opera- 
tors 

F b 1 ,...,b m ;a 1 ,...,a n ^ ■ ■ ■ ^ m \9l, • • • ! #n) 

= bl ,..., V (9[, . . . , ?J*(O)|0i, . . . , 0„) oi ,...,a„, (51) 

or form- factors (FFs), can be computed exactly 47 . The 
related perturbative approach is called Form- Factor per- 
turbation theory (FFPT)2£. 

According to Ref. [2^, the first order term in the ex- 
pansion of the mass variation of the particle A a is given 
by 

5m 2 a ~ 2g „<0|tt(O)|0> a = 2g C^' F^ a (m) , (52) 

where the two-particle FF, F* Oa (0i — 9 2 ) = 
ai,a 2 depends on the difference of the 
rapidities only if \1/ is a scalar operator. In (|52ll we used 
the crossing property valid for Lorentz scalars 

■ F 6*...,6 m ;oi,...,o„(^i) ■ ■ -jflmJ^l) ■ • • > ^n) (53) 
m 

= Y\_ Cbl ' bl F V 1 ,...,V m ,a 1 ,...,a„0'li ■•■^m^l!-'. On) 
J'=l 

where 9 = 9 + m and C is the charge conjugation ma- 
trix satisfying two requirements: C l — C, C 2 — I. For 
Majorana fermions C is trivial because they are neutral. 

Since this is a strong coupling (IR) analysis, if g in 
(|5U|) scales under RG, it has to be replaced in l(5^|l by its 
renormalized value at energy of the order of the largest 
mass in the theory 

g ^ g eff~ g (m). (54) 

The definition of the coupling constant at energy m is 
somewhat ambiguous. The problem is that RG equations 
are universal only in the first loop and beyond this they 
depend on the regularization scheme. By introducing 
g e ff we stretch these equations to their limit. Therefore 
we will not be able to establish a rigorous relationship 
between the IR and UV parameters. 

In analogy with QM degenerate perturbation the- 
ory, the perturbed masses for degenerate multiplcts 
are obtained by diagonalizing the matrix {M m „} = 
{F* m (i7r)}, where indices n, m belong to the degenerate 
multiplet. If the symmetry of the perturbing operator is 
less than the symmetry of the multiplet, the perturbation 
will split it. 

Assuming that the IR coupling constants in l|52|l are 
smooth functions of the bare ones and F*(i7r)o, is finite, 
the spectrum evolves adiabatically. However, there are 
effects which do not appear in the first order. The pertur- 
bation may also generate an effective attraction between 
the particles, and lead to an enlargement of the particle 
content of the theory through the formation of bound 
states below the two-particle threshold. For small g e '* 
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the appearance of bound states can be studied solving 
the eigenvalue equation for the two-particle wave func- 
tion. This will be done in some detail later in the text 
(cf. ©). 

A more dramatic change in the particle content of the 
theory happens when if? is non-local with respect to the 
operator, ip a , that generates the particle o2£. In this case, 
as shown below, the two-particle FF, F^ a (0), has a pole 
for 9 = in and then the mass variation l|52[l diverges. This 
infinite mass variation implies that the original particle 
a (it will be called a "quark" in analogy with high energy 
physics) is confined and disappears from the spectrum of 
the perturbed theory. The resulting spectrum consists 
of quark bound states - "mesons" . Their masses may be 
higher than the two quark threshold. The analysis of the 
meson spectrum will be done in Sec. IVII Dl 

The way the non-locality properties of an operator af- 
fect its FFs can be seen as follows. From the general 
theory of FFs it is known that the n-particle FF has 
kinematical poles, associated to annihilation processes, 
whenever the rapidities of a particle and anti-particle dif- 
fer by in. The residue is given by^£ 

• A)= (55) 



1 



a 27rl 7a 




where 7 a .^ is the semi- locality index between ip a and ^f>. 
It is easy to see that for two particles (n = 0) the right 
hand side vanishes if 7 Qi * is an integer. Then a two 
particle FF has a pole at 9 = in only if the operator is 
non-local. As a consequence of l|55[) the mass variation 
(|52[l induced by a non-local operator is infinite. From 
(|55l it follows that for (9 ~ 9') 



5a,a'5b,b' , 5 a ,b'5b,a' a 

—^b a ,b{2,9 ) 



-9') 2 {9 + 9'Y 



, (56) 



with a a . b = (1 - c 2l7r ^.*)(l - e 21 ^".*). This will be 
needed in the following. 

The factorized scattering approach can also be con- 
structed for massless integrable models^, despite the 
subtleties in defining a scattering theory between mass- 
less particles in (1 + 1) dimensions. In this case the 
excitations are right and left moving particles, Ar(6), 
Al(9), with dispersion relation en(9) = pr{9) = 4fe e 
and ez,(0) = ~pl(9) — ^[e~ e , where M here is just a 
scale (for simplicity we consider only particles with no 
internal indeces, this is the only situation that we will 
encounter in the following). Also the FFs can be defined 
and computed in analogy with the massive case^i, and 
FFPT can be used to study the mass generation that, to 
first order, will be given by^ 

5m ~ g F$ L (in- oo), (57) 



where F% L {9) is the right-left FF 



|¥(O)|4(0i)4(0 2 )>. 



(58) 



Also in this case confinement is related to the non-locality 
properties of the perturbing operator—. 



B. Non- locality properties of the perturbing 
operators 

In order to apply the above methods to deformations 
of GN models like in I(25I28|) . let us study the locality 
properties of the perturbing operators with respect to 
the GN particles. 

Consider first a massive O(N) GN model (N > 2). 
Generally speaking an operator of the form vf^ = V'jjV'x 
is local with respect to the vector particles, Af i; but not 
with respect to the kinks, A^. This can be easily seen 
even without explicitly introducing kink creating opera- 
tors. In fact the kink is an elementary excitation of the 
GN model that interpolates between positive and nega- 
tive minima of tpR^I ano - therefore this operator changes 
sign across a kink configuration. This implies that the 
non-locality index between the kink creating operator 
and is e 21 ^ 1 — — 1, but, at the same time, an operator 
that is a product of two fermionic bilinears belonging to 
the same GN model, vE^ = ^r^l^r^li is local also with 
respect to the kinks. As a consequence, Eq. l|52"|) can be 
safely used to evaluate the first order effect of deforming 
the 0{N) GN according to (HJ, since the perturbation 
is local with respect to both kinks and fermions. This 
implies that the spectrum will evolve adiabatically, the 
only important effect is to split the degeneracy between 
some states as discussed in the next section. 

On the other hand, if one takes two different massive 
GN models, say # (;v)[x] and H 0(M) [£] (N,M > 2) 
and couples them like in 1)28(1 . the perturbation, ^3 = 
XrXl£r£li i s non-local with respect to the kinks of the 
two GN models, thus leading to their confinement. This 
is in agreement with the semiclassical analysis of Sec.lVI 



C. First order effects: Perturbations around the 
0(6)-symmetric point 

We first apply FFPT to study the spectrum of vec- 
tor particles around the 0(6) symmetric point. We con- 
sider small deviations of coupling constants in from 
the point where all of them are equal g a = g + 5g a 
(a = [p, — ], [c, ss], [a, —], [c, at], [<r, +]). The perturbation 
is local and then Eq. I|52|) can be safely applied, with 
g a replaced by g a (m). It should be emphasized that the 
coupling constants g a carry information about physical 
quantities only when they are small. Beyond the first 
loop approximation RG equations are not universal de- 
pending on the regularization scheme. We assume that 
the interaction between physical particles do not change 
substantially once g a ~ 1 and therefore stretching the 
RG equations to the limit of their validity ln(e/m) ~ 1, 
we can estimate the anisotropy of the coupling constants 
in the strong coupling regime. This assumption is not 
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justified but its consistency with know results will be 
checked. 

For Sg a /g <C 1 we can linearize the RG equations 0. 
The details of the calculations with all notations are given 
in Appendix IdI We find 

Sg a (e) = A [log(e/m)]- 2 + (59) 
(T a , 1 C_ 1/2 , + + T a , 2 C7_ 1/2 ,_) [log( £ /m)] 1/2 

+ (r„, 3 c 1/2j+ + r , 4 c 1/2 ,_) [iog( £ /™)]- 1/2 

where C g ,± ~ S q g(0)gQ q (S q stands for a particular linear 
combination of couplings) are RG invariants, that can 
be obtained explicitly using the results of Appendix IDI 
and A ~ ^Sg a (0)/gQ. At scales of the order of the 
mass we can set ln(e/m) w 1 in (|50|> and substitute the 
result into Eq. I|52l) . By dimensionality considerations the 
matrix element of the current-current product is const 
m 2 . We remark that this normalization yields the correct 
estimate of the mass change at the uniform variation of 
the coupling constants. Indeed, in this case, according to 
(|59|l , the change of the effective coupling at energy m is 
Sgo/g 2 , and ljS*2"|l gives 8m ~ Sgo/g^m. On the other hand 
we obtain the same estimate from the known dependence 
of the mass on the bare coupling: 



2 dm 2 2m 2 n 
dm = — — 8 g = „— 8g 



dgo 



9l 



(60) 



It is interesting to estimate the mass splittings for the 
case of small 50- Then the largest RG invariants that 
produce a splitting of the multiplet are Ci/ 2 (the most 
relevant term, proportional to Aq, is the same for any 
fermion). Taking only them into account we obtain the 
following mass corrections: 



)m c 



-) 



8m 



m 

Cl/2, 

0) 



[8g p -(m) + 8g c 



,(m) + 38g CtSt (m)} 



rn 

~ C1/2, 

8mf r ^ 



[38g a -(m) 
- 2C!/ 2: _ 



28g c ,ss(m)] 



m 



[28g a , + (m) + 8g a -(m) + 2dg c>8t (m)] 



"Cl/2H 



(61) 



One can conclude from Ea. (|61|l that for any finite RG 
invariants the 0(6) vector multiplet is split. Nevertheless 
the splitting is very small in the regime of validity of this 
analysis. 

Let us apply now Eq. I|52|) to study the evolution of 
the kink masses. Since they are degenerate we need to 
diagonalize the matrix with elements ki (Oi | ^r^l^r^l I 
Oj)kji where fcjj indicate kinks of the degenerate multi- 
plet. The needed FF has the form 1 



52 



F, 



a . b 



(%) = (0 1 



0—8 a<b )C ki<k J{6), 



(62) 



two-kinks continuum 



2(M,+ M 2 ) 



O(M) kinks M 2 

, fermions 
0(N){ 




M J+ M 2 H 



kinks M 



bound states of O(M) kinks 



bound states of O(N) 
kinks 



adiabatic evolution 
of fermions 



0(N)xO(M) 



FIG. 2: A qualitative picture of the spectrum of the O(N) x 
0{M) GN model. 



where C is the charge conjugation matrix and the func- 
tion f{6) is finite at 9 = m. From this, using (|53|) . one 
gets that, for b, 

uM I ^1^1 I 0j) kj = S ki , k .f(9), (63) 

and then the matrix is diagonal. Using now Eq. I|52l) one 
finds that the mass variation is the same for any kink and 
then the multiplet is not split. 



D. Beyond the first order 

Let us now turn to the problem of kink confinement 
induced by a non-local operator. This analysis in im- 
portant to study the spectrum close to the first order 
transition line. In this case the FF in (|52|l is infinite 
and this formula cannot applied. Nevertheless, for g e ^ 
sufficiently small the lower part of the meson spectrum 
can be studied within the two-quarks approximation de- 
scribed in Ref. |48| in the context of the Ising model. It 
is convenient to parameterize the ZF operators directly 
in terms of momentum; the relationship between the two 
parameterizations can be found in (|48|l . The two-particle 
wave function can be written in the center of mass frame 



-Hoc. 



d^ b { P ) A{(p)Al(-p)\0). (64) 



In general a and b can be different particles with masses 
M a and Mf,, but both non-local with respect to ^. Please 
note that the wave-function is not written in an explicit 
relativistic invariant form although the states have a rela- 
tivistic normalization. The eigenvalue problem gives the 
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following equation 

[E - e a (p) - e b (p)} ('ljja : bip)Sa,a'Sb,b' + lpb,a(p)Sa,b>Sb,a>) 
= 9 £ff / dp' a , b ( P , -p | *(0) | p', -p')a>,b> 4>a',b> tfffiS) 



where g e ^ is the renormalized coupling constant at en- 
ergy ~ (M a + M b ). For energies close to the threshold 
(small momentum) one can approximate E — e a (p) — 
e b (p) ~ e - p 2 /2fi, with e = E - (M a + M b ) and 
fi = M a M b / (M a + M b ). If the operator is non-local, it 
follows from (|55|l that, the right-hand side is dominated 
by the double pole in the four-particle FF. It will be con- 
venient for us to distinguish between two cases. One is 
the confinement of identical particles and the other is 
when confining particles are not identical having, in gen- 
eral, different masses. These two cases are qualitatively 
similar, but there are differences in details. 

We consider first the confinement of two kinks that 
belong to different GN models induced by the operator 
^3 = XrXl^r^l- Let us call ai kinks that belong to 
Ho(n)[x] an d bi kinks that belong to Ho(m){£]- Given 
the structure of the potential, the FF factorizes as 

ai,6i(Pl,P2|*3|p / l.P / 2)o2,6 2 ( 66 ) 
=«i (Pl\(X R X a L)\p'l)a 2 MtfR&M)*. 

According to (|55|l and (|56|l for p ~ p' we can approximate 



aAPiliXRXL^p'i) 



2iM a 



(Pi - Pi) 



jt{XrXl)S 

ai ,a,2 

(67) 



where we have used the fact that e 2l7r7a '*3 = — 1. Clearly 
the same approximation is valid for b 1 {p2\(^R^L)\P2)b 2 - 
From this we get the eigenvalue equation 

/OO 1 / 

-, _ p i>aAp')m 
-co IP P ) 

were V indicates the principal part and 

Kb = K f biXRXL){^L)M a M b . (69) 

After a Fourier transform to configuration space Eq. I|68|) 
takes the form 



(-e+A , b |X|- 



1 d 2 
2n dX 2 



)i> a , b (X) = (70) 



where 



X) : : / ^(> X V„, 6 ( P ). 



(71) 



In general this equation must be supplemented by sym- 
metry conditions. 

The eigenvalue equation (|71|l describes a particle con- 
fined in a linear potential. The solution is well known 49 . 



In the region X > one can introduce the new variable 
C = (2A a , 6 (J,) 1/3 (X - e/X a .b) so as to yield 



d( 2 



0*(0 = o 



(72) 



at £ > — (2^i/A 2 & ) 1 / 3 e, from which it follows that the 
solution is the Airy function 



^ Qjfe (C)cxAi(C). 



(73) 



The eigenstates are determined either by the zeros of the 
derivative of the Airy function Ai'(£) at X = (symmet- 
ric wave functions), or by the zeroes of the function itself 
(antisymmetric wave functions). Calling them Q we have 



E t - (M a + M b ) = (\l b /n) y % ; 
(Ai'(-Ci)=0 or Ai(-Ci)=0). 



(74) 



Clearly, if the particles a, b are identical the above anal- 
ysis requires minor modifications as we will see below. 

In agreement with the semiclassical arguments of Sec- 
tion IV we see that the kinks of the 0(N±) and O(A^) 
models disappear from the spectrum and are replaced by 
the generations of mesons which energies belong to the 
original two-particle continuum starting from Mi + M 2 . 
Stable meson states must have spectral gaps smaller than 
Ei < 2Ei w 2(Mi + M 2 ). The above analysis is quantita- 
tively valid if the number of meson generations is large. 
For the simplified model 0(3) x 0(3) it is easy to ex- 
press this number in terms of the RG invariants. Setting 
Ai.2 =constMi.2 we find from Eqs. H89|l below that at 
g + ~ — 1 the effective coupling constant is 



[2|C 1/2 |+C_ 1/2 ]- 2 



(75) 



provided \Cx/ 2 \ + C-1/2 >> 1- Then using Eq. JHSJl we 
obtain 



G?y/i) i/3 ~(M 1+ M 2 ) 



Mi + M 2 



2/3 



(76) 



Since Q 



,•2/3 



at i >> 1, this sets the limit for the 



number of meson generations as 



9eff 



} f (y/Mi/M 2 + ^M^/Mx 



(77) 



Thus it appears that the limit of weak confinement 
can be achieved by either taking small ff e // or consid- 
ering quarks with vastly different masses. These two 
routes however lead to noncquivalent limits. The limit 
M1/M2 — * is different from the limit g e ff — > and 
cannot be studied using the above formulae. The reason 
is that when the confinement radius (g e //AfiM 2 /i) -1 / 3 
becomes of order of the size of the lightest kink M-j~ , 
one can no longer use Ea. (|67|l . Thus the above analysis 
is valid only at 



Mi » g e f/M 2 



(78) 
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At smaller mass ratios one can either use more accu- 
rate expressions for the FFs or resort to the semiclassical 
analysis of Section TV Bl As we know from this analysis, 
confinement of a light particle on a heavy kink produces 
a zero energy bound state and the scattering continuum 
separated from the bound state by a gap of the order of 
g e ffM 2 . This is of the order of the energy of the lowest 
bound state (gl b /fi) 1/3 at M x ~ g e ffM 2 . 

Let us repeat the above analysis for confinement of 
identical particles (M a = M b = M) induced, for instance, 
by an operator 'P. As we discussed in Section IV, this 
leads to formation of vector particles. Within the integral 
of lt)5|) . we can approximate 



a , b {p, -p\ *(0) \p',-p')a',V 



{p-p') 2 (p + p'f 



(79) 



With few manipulations Eq. (|fi5[) can be put in the form 
JZOl with A a , fc = 4g e ^M 2 (^). In this case the particles 
are indistinguishable so only antisymmetric Airy func- 
tions are acceptable. The summary of the above results 
is shown schematically in Fig. 

It should be also noticed that if the unperturbed the- 
ory, like in the case of the O(N) GN model, possesses 
bound states of mass M a j, between particles a and b (as- 
sociated with poles of the S-matrix on the physical strip), 
the eigenvalue equation Ij65(l has to be modified taking 
into account this residual attraction 

[E - ei(0) - e 2 (0)] -+[E- ei(0) - e 2 (0) + AE a , b ] (80) 



where AE a ^ b = Mi + M 2 
Eq. l(70"|) is modified as 



M a &. As a consequence, 



1 d" 
2fi dX 2 



+ \\X\-bo5(X))il> a ,f,(X) = 0. (81) 



The delta function is chosen in such a way that vector 
particles appear also at A = 0. For the 0(3) case, &o = 0. 

The formation of bound-states induced by a local per- 
turbation can be discussed in similar terms. For small 
g e ^ the two-particle wave function can be taken in the 
form l|64|) from which the eigenvalue equation I|tj5fl fol- 
lows. The crucial difference with respect to the case 
discussed above is that, since the operator is local, the 
4-particle FF in l|(j5|) only have first order poles giving 
rise to a non-singular potential. Bound states form if the 
perturbation is such that there is a solution of (|rJ5|) with 
energies lower then M\+M 2 . This is what is found, for in- 
stance, for the 0(4) GN perturbed by the 1 2 anisotropy, 
as discussed below. 



E. Coupling of massive and massless modes 

The situation close to the U(l) and 1 2 critical points, 
described by Eq. (|28[1 with N 2 = 1, 2 is more tricky since 



some gapless modes are present. This problem was al- 
ready addressed in Sec. IV Gl and FFPT cannot add much 
to that analysis. Nevertheless we find it instructive to 
briefly discuss the problem also in this framework. Let 
us look for simplicity at (|28|l with N\ = 3, N 2 = 1 and 
uniform deformation 5g ab = 5g. For Sg = the spectrum 
consists of a massless Majorana mode and massive 0(3) 
kinks. According to l|52|) and l|57|l the mass variation is 
given by 



CO 



Sm 3 ~2Sg (xW)^'(i 
6m 2 a ~2 6g(e R eL) a (e\x a R xl\0)a 



(82) 
(83) 



where the indices 3 and a indicate massless Majorana 
fcrmions and massive 0(3) GN kinks respectively. The 

VEV (XrXl) i s different from zero and F R ^ ' (9) is con- 
stant for any value of 8, then from (|82|l it follows that 
the perturbation induces a mass linear in Sg, i.e. as soon 
as the coupling is turned on a massive singlet appears in 
the spectrum. 

One needs to be careful in the analysis of the second 
equation. In fact vanishes in the unperturbed the- 

ory, implying that the mass variation of the kinks van- 
ishes in the first order. Nevertheless (£fj<!;£) ~ Smz and 
then (|83|l can still be applied, just keeping in mind that 
the effect on the kinks mass is of higher order in 8g. 
Since the operator (xrXl) * s non-local with respect to 
the kinks, the rhs of Eq. (|83|l diverges and kinks con- 
fine. One should notice that besides the one considered 
in H83f) . there are other second order contributions to the 
mass variation of the kinks, nevertheless they will not be 
able to compensate the divergence and will not modify 
qualitatively the kink confinement. It can be shown that 
the confined states will always be unstable and then, in 
agreement with the semiclassical analysis, the only stable 
particles will be singlet and kinks with quantum numbers 
of the 0(N + 1) model. 

These results can also be checked using the opposite 
limit Sg ~ g\ = g 2 and treating the model (|28|l as a 
perturbation of the 0(4) GN (remember that for N = 4 
the GN model has only kinks in the spectrum). This 
perturbation is local with respect to the 0(4) kinks and 
then produces an adiabatic change of their masses. At 
the same time we know that the 0(4) model also contain 
(unstable) vector particles above the two-kink threshold. 
As described in the previous section we can use FFPT 
to check whether the perturbation induces further at- 
traction that stabilizes them. Since the vector particles 
are bound states of kinks, for small deviations from the 
0(4) point we can write their wave function like 16411 . 
with a and b being 0(4) kinks, and study the stability 
of the bound state by solving the eigenvalue equation 
(|65|l . Even without entering the details of the compu- 
tation one can easily see that the form of the potential 
naturally split singlet and triplet bound states and in par- 
ticular the interaction between kinks in the singlet state 
is attracting while the one between kinks in the triplet 
state repulsive. Then it turns out that the interaction 
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2-kinks continuing 



O(N) fermions 



O(N) kinks 



0(N+1) 
vector multiplet 




0(N)XZ 2 



0(N+1) 



1) model presented in App. [U] The spectrum consisting 
of 5+1 vector particle multiplet and a quartet of kinks. 
The symmetry of the spectrum remains 0(5) XZ2, which 
corresponds to the splitting of the vector multiplet as 5 
+ 1. Following the arguments of Section IVlI CI we find 
that at large C >> 1 the mass ratio of the different 
vector particles is ~ gx(m) ~ C~ 2 / 3 . At C << 1 the 
vector multiplet approaches the 0(6) symmetry with the 
splitting ~ C. 

Then in general, the SCd phase has kink particles with 
the 0(6) quantum numbers and vector multiplet con- 
taining a triplet, a doublet and a singlet with different 
masses. All masses have the same sign. By going trough 
Z 2 QCP the singlet mass changes sign and one finds the 
Orbital Antiferromagnet (CDWd) phase. 



FIG. 3: A qualitative picture of the spectrum of the model 
interpolating between the O(N) x Z 2 and 0(N + 1). GN 
model. 

stabilizes the singlet that will appear in the spectrum. 
This implies that the singlet state will be stable for any 
value of Sg < g 3 and will cross the threshold at the 0(4) 
point (see Fig. 0. We can repeat the same procedure 
for Ni = 5, the main difference is that the unperturbed 
theory has fermionic states in the spectrum. 



VIII. THE SPECTRUM AROUND PHASE 
BOUNDARIES 

The analysis of the previous Sections allows us to un- 
derstand qualitatively, but rigorously, the spectrum of 
the model © close to the phase boundaries, where it 
can be seen as a perturbation of some integrable models. 

A. From SCd to CDWd (SCs to CDW) phases 
through Z 2 QCP 

At the Z2 QCP the Hamiltonian J5J is described by 
a deformed 0(5) GN model decoupled from the Critical 
Ising model. The effect of the anisotropy is to lift the 
degeneracy between the vector particles. Away from the 
Z2 QCP the coupling is 

- 2(&£!) {9,-(XRXl) + ffc„[«ki) + ■ (84) 

It is convenient to study the spectrum for g a + = g p ._ = 
9c,st = 9; 9a- = 9c,ss = 9x when the 0(5) symmetry is 
exact. The RG equations 

9 = -3g 2 -g x , 9x = -^9xg (85) 

have the following RG invariant: C = (g 2 — g^/g^ 2 . 
To understand the spectrum we can use results from the 
previous sections as well as the exact solution of the (5 + 



B. From one SC (CDW) phase to another one 
trough U(l) QCP 

At the U(l) transition two Majorana fermions are 
massless and the rest is a deformed 0(3) GN model 
equivalent, for small deformations, to the (3+l)-modcl 
discussed in Sec. IVII Thus the spectrum consists of the 
Gaussian massless mode, the SU(2) kink doublet and a 

singlet particle of a lower mass. The Q { ~ ] field acquires 
a mass gap as soon as the coupling between the two sec- 
tors is on. The effective field theory describing the lowest 
excitations is the sine-Gordon (SG) one: 

s = i[0 M e£->] a + «oosb9e£->], (86) 

with v =_i([<7 SlCt (xM + 9c,ssiUl]) and 1 = 4tt(1 - 
<7 Pi _/27r) . The SG kinks with period 2ir/f3 constitute 
a low-lying doublet. Thus the vector multiplet is re- 
duced to the U(l) doublet and the Majorana singlet. For 
(3 2 > 4-7T (repulsive interactions) there are no sine-Gordon 
bound states, but if (3 2 < An, are also additional bound 
states (breathers). These particles are absent in the 0(6) 
GN model. As far as the kinks are concerned, they ac- 
quire both spin and U(l) quantum numbers via creation 

of bound states between half (anti)kinks of <di ' with 
period it/ (3 and the SU(2) kinks. 



C. Around the first order line 

At the first order transition the spectrum consists of 
kinks with different masses, Mi and M2 belonging to 
the isotropic and anisotropic (£7(1) X Z2) 0(3) GN mod- 
els respectively. As soon as the interaction is turned on 
all kinks confine and one has fermions close to the two 
kink thresholds 2M\ and 2M2 and new kinks close to the 
M1+M2 threshold. As we discussed earlier in Sec. IVII Al 
the region of weak confinement where one may observe 
several generations of kinks, is likely to be narrow. 
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FIG. 4: RG flows JHHl for C x /i 
the flow with C_i/2 = 0. 



0. The thick line represents 



FIG. 5: RG flows lHHJ for C 1/2 > 0. 
the flow with C_i/j = 0. 



The thick line represents 



IX. FURTHER FROM THE BOUNDARIES: 
0(3)xO(3) MODEL AS A SIMPLIFIED CASE. 

With the results of the previous sections at hand we 
will return to the case N = 6 and consider the model @ 
with a more symmetric form of the interaction 



V = -g + (x a RX a L) 2 
-2gx(x a R xm b R e L ) - g- 



(87) 



as it was done in lUloSlsd Though the subsequent RG 
analysis has a significant overlap with the one conducted 
in these papers, we include it here for completeness. The 
RG equations have the following form: 



9+ = 9% + Zdxi 9- 

gx = 2g x (g++g-), 



(88) 



In this case we have managed to find the RG invariants 
cxplicitely 



C_i 



/2 



5+5- 



V\9X\ 



/2 



5- 



2V\9x\ 



(89) 



For g a = g + 5g a these invariants correspond to the ones 
found in Section IvTl"C1 If the 0(3) x 0(3) model is re- 
placed by 0(N)xO(N) the power 1/2 in JHU is replaced 
by (N-2)/(N- 1). This means that in the large TV limit 
the ratios of coupling constants become RG invariants 
which coincides with conclusions of Section We also 
observe that a sign of the quadratic form 7 is RG invari- 
ant which gives support to the 1/N analysis of Section 
IVl Using these RG invariants one can integrate Eas. l|88|l . 
It should be kept in mind however, that equations l)88|) 
together with the explicit form l|89[) are valid only at 
\g a \ < 1. As soon as one of the couplings becomes ~ 1 
the scaling of this coupling should be stopped and the 
problem reconsidered. 



The solution of this equation is given by 
rax(t) dg 



At = ± 



<w(0) g^C^ 1/2 g l / 2 + Cl /2 g + g 2 ' 



(90) 



that can be expressed in terms of elliptic functions of 
complex modulus and argument. 

The phase diagram has three areas corresponding to 
different properties of the interaction tensor. 

• The signature of the interaction tensor is positive: 
C1/2 > and C_i/ 2 > 0. All couplings scale to 
zero. 

• The signature is negative (Ci/ 2 < 0) and C_i/ 2 > 
0. The system scales to strong coupling. The en- 
ergy scale on which the spectral gaps are formed is 
given by 



M ~ A exp 



dg 



l«(0)| 5 v /C-i/25 1/2 + C 1 2 /2 5 + 5 2 



;9i) 



At C_i/2 > the integral in the exponent diverges 
at gx(0) — > 0. Therefore the theory has a proper 
scaling limit M =const gx(0) ~~ * 0,A — > 00. This 
limit is characterized by the RG invariants C±i/ 2 - 
The case C±x/2 = corresponds to the 0(6) GN 
model. The opposite limit max(C-|- 1 /2) >> 1 de- 
scribes the 0(3)xO(3) model in the regime of weak 
confinement (see Sections IV BIVIID|) . 

C_i/2 < 0. The system scales to strong coupling 
for any sign of the bare diagonal couplings. In this 
case the scaling trajectory for gx (t) bounces off its 
minimal value g* given by the root of the equation 



? /2 + C 2 1/2 [g* 



1/2 



\C-i 



/2I 



(92) 
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If 9x(0) > the weakest interaction is achieved at 
the energy scale 



m a /m 



E* = A exp 



9x(0) 



dg 



g^C- 1/2 g 1/2 + C 2 /2 g + g 2 

and the strong coupling regime is reached at 



93) 



M = E* exp 



dg 



a' g^C-i/ 2 g 1/2 + C 2 1/2 g + g 2 _ 



(94) 

As we said in Introduction, the maximal spectral 
gap corresponding to gx(Q) —> +oo is 



M ~ E* 2 /A 



(95) 



Since integral (|§U)l never diverges at C_xli < 0, 
the RG time in which the strong coupling limit is 
reached always remains finite. Therefore at finite 
C_xj2 < the scaling limit does not exist. 

If we relax conditions on the scaling limit we can ven- 
ture into the area C_i/ 2 < 0, as far as the spectral gaps 
remain much smaller than the bandwidth. For this we 
need g* to be small which is achieved by making C\ small. 
Thus the effective low energy theory describing the state 
with emergent attraction is the C-1/2 = field theory 
with non-relativistic corrections in M/A. The condition 
C_i/2 = does not put any restrictions on C1/2; m the 
realistic two-leg ladder, where bare lattice interactions 
are not small, C1/2 may have any value. Ea. (|89|l point- 
ing to small values of C±\/ 2 at small bare couplings is 
not valid in that limit. 

Taking the above into account we can qualitatively de- 
scribe the overall spectrum of the 0(3)xO(3) using re- 
sults of Sections IVII CIVII Dl The results are summa- 
rized on Fig. 



A. The mass is induced in the otherwise massless 
sector 

It remains to study the situation when some coupling 
constants are much larger than the others such that they 
reach the strong coupling limit first. This analysis is 
complementary to the one of Sees IV CI and IV11 hi In 
this case one should study the spectrum using two-stage 
RG. There are two cases to consider: 

• One of the sectors is massless at gx = 0. Let it be 
the + sector: #-(0) < 0,.g + > 0, so that at gx = 
the interaction in Xa sector is marginally irrelevant. 



The condition C_ 



1/2 



< is fulfilled. As we shall 



demonstrate, the vector particles do appear in this 
sector at gx 7^ 0, but the symmetry is not restored. 

• Both sectors are massive at gx — 0. Let < 

\gx(0)\ « g+(0)g-(Q) and -g+(Q) ~ - 5 _(0) > 0. 




*- 1/C 



FIG. 6: The qualitative picture of the mass spectrum of the 
0(3)xO(3) GN model as a function I/C1/2. The dashed lines 
show masses of mesons. 



Then the coupling gx{M) is still weak. Following 
the arguments of the previous sections, we conclude 
that the interaction leads to the creation of gener- 
ations of 0(6) kinks with masses close to Mi + M 2 
and formation of two generations of the vector par- 
ticles: one with masses close to 1M\ and the other 
one with masses close to 2M 2 - 

We use the two stage RG. Let M3 is the mass scale of 
the 0(3) GN model in the — sector. At this energy the 
renormalization of g^ coupling terminates and a nonzero 
average is formed: 



i(d(^ 3 _1 )ei(o)) « -m 3 



(96) 



Replacing operator by its average in the low energy 
effective action for \ fermions we obtain the effective ac- 
tion like 1331) with the UV cut-off M 3 : 

^e// = ^XaJfidftXa + iM XaXa ~ g+(M 3 )(XaXa) 2 (97) 

where M = i5 X (M 3 )(^(M 3 " 1 )^(0)). The form of this 
Hamiltonian suggests that the vector triplet in the — sec- 
tor is stabilized. For the procedure to be self-consistent 
the mass of this vector particle should be M + « M3. 
Let us assume that g + (M 3 ) is still positive, which is the 
case at sufficiently small gx(0)- Then the downward 
renormalization of g + will continue till the energy M + 
(the resulting mass for x fermions). The latter is deter- 
mined by the equation 



M+ = 



M 



1 



g+(M- L J 

47T 



ln(M 3 /M+) 



(98) 



We see that the presence of the repulsion in the + chan- 
nel can further reduce the mass of the vector particle. 
The model (|97f) was introduced in Ref. |^ to describe 
low-energy properties of the spin S=l Heisenberg anti- 
ferromagnetic chain. 
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B. 4 +2 splitting 

The trickiest and at the same time the most realistic 
case is when the modes split as 2+4 with repulsive in- 
teraction between the four magnetic modes. In this case 
neither group can produce a mass gap and the gap is 
generated solely by the interaction. This happens, for in- 
stance, in the model with no inter-chain tunneling in the 
presence of the density-density interaction. Then only 
the 2kp components of the density interact leading to 
the arrangement 

9a,+ = 9a,- = v>0, g CiSS = g c<st = g 

where v originates from the intra-chain exchange inter- 
action. The 0(6) symmetry is split as 0(4)xU(l) 



igcospei-^XRXl+v^XRXi) 2 - 



(99) 



The RG equations 



9 P - 



-4g 2 , v = -2v 2 -2g\g= -g{g p ,- + 3v) (100) 



where /3 2 /4tt = l-\-g Pi _/2-K. In the realistic case of repul- 
sion g p _ < the system scales to strong coupling with v 
changing sign. Assuming that —g p — is the largest among 
the coupling constants and neglecting its renormalization 
we arrive to the approximate solution 

g = g(Q)e^-\*, v = v(0) - f^l( e 2 l^-l' - 1),(101) 

\9p,-\ 



such that v changes sign at 



1 



hi 



\9p,-H0) 



3 2 (0) 



(102) 



which fixes the upper limit for the gaps E* 2 /A ~ Ae -2 * . 
This corresponds to the scenario when the mass is gen- 
erated by the effective attraction among the fermions. 

At small d — (3 2 /Att < 1/2 another scenario can be 
realized. To understand the spectrum we resort to 1/N 
approximation where N is the number of fermion species. 
Assuming that the fermion masses are larger than the 
boson ones, we integrate over fermions and obtain the 
following effective potential for Oc = <f> held: 



V 



2 JV 



4tt 



cos 2 (/3^>) In 



A 



S|cos(/30)| 



(103) 



which can be approximated as a pure cosine. This po- 
tential is relevant only when d < 1/2 and for d > 1/4 the 
spectrum contains only kinks. Let us call their mass gap 
77i0. Then the estimate for the fermion gap is 



m x w ff (cos(/30)) = gC l3 k{m <t ,/K) d 



(104) 



where Cp is a constant. The fermion gap is indeed much 
larger than the kinks gap which justifies the integration 



over the fermionic modes. Thus the spectrum in this case 
consists only of the 0(6) kinks. The kinks acquire quan- 
tum numbers from the spinor representation of the 0(6) 
group through attachment of the fermions zero modes 
(see Section W} . 



X. CORRELATION FUNCTIONS AND 
EXPERIMENTAL PROBES 

A detailed analysis of correlation functions is outside 
the scope of this paper, so we restrict this discussion to 
some qualitative remarks. 

Among the available experimental methods the ones 
which probe correlations functions at various frequen- 
cies and momenta are Angle Resolved Photoemission 
(ARPES), inelastic neutron and X-ray scattering. All 
other probes measure either local correlation functions 
(Nuclear Magnetic Resonance (NMR), tunneling spec- 
troscopy) or correlation functions at zero momenta (op- 
tical conductivity, Raman scattering). 

The only known system which effective Hamiltonian 
resembles the one for a single two-leg ladder are single 
wall carbon nanotubes. However, though the doped nan- 
otube is described by the effective Hamiltonian (0!—, the 
values of the bare couplings for the nanotubes of available 
size are such that the gaps are extremely small. In other 
experimentally available systems, such as the telephone 
number compound mentioned in Introduction, ladders 
are packed into a three-dimensional (3D) arrangement. 
In that case interactions of gapless charge modes from 
different ladders may lead to 3D ordering. 

Let us discuss the behavior of correlation functions 
above 3D phase transition (if such transition occurs). In 

this case both <&c + ^ and 9^ + ^ exponents have power law 
correlation functions. Since these exponents enter the 
majority of experimentally measurable correlation func- 
tions, emission of all massive excitations are accompa- 
nied by emission of gapless bosonic modes. This makes 
the massive particles incoherent. The best one can do at 
the circumstances is to look for operators whose correla- 
tion functions contain matrix elements corresponding to 
emission of just one massive particle (and, of course, a 
cascade of gapless bosons). 

Having this in mind, let us consider, for example, the 
SCd phase. In this phase all a (/x) fields are locked 

together with cosf^/TrOi ] (sin). Then the operators 
with matrix elements between the vacuum and a single- 
particle state of the doublet vector particle are the 
Fourier components of the particle density at 2k^ : 



P {2k^) = Y.Kv L ^ 



(105) 



ie v c e v <= [o-i 020-3 cr T M1M2M3M0J 



As far as the magnetic triplet is concerned, it is emitted 
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at 2(k ( p> + k ( p>) by the magnetization operator 

S a (2k F ) = R+^a^L^ = (106) 

zira 

where 

N a = (AilM20-3,Ml°'2Ai3,Cri/i2Ai3), 
M a = ((Ti(T2M3,0'lM20'3-MlO-20'3) 

The corresponding quasi-coherent peaks disappear below 
T c since the g£ + ^ is locked and correlations of e L ^® ( <= +) 
become short-ranged. However, if the SCd phase is 
replaced at stronger interactions by the Akp ordered 
Wigner crystal (recall the discussion in Section II B), it is 

$>c held which is locked. Then the all the above peaks 
become sharp (this was in fact observed in the telephone 
number compound^). 

XI. CONCLUSIONS 

In this paper we made an attempt to outline the picture 
of the excitation spectrum of a field theory with a sym- 
metry of a doped two-leg ladder, namely the anisotropic 
0(6) Gross-Neveu model. This model is not integrable 
and then it is not possible to obtain exact results ex- 
cept for some specific points in parameter space. We 
combined information coming from the RG analysis with 
semiclassical methods and form-factor perturbation the- 
ory. It follows from our analysis that throughout most 
of the phase diagram the spectrum consists of degener- 
ate quartets of kinks and antikinks and the multiplet of 
vector particles split as 3 + 2 +1. This basic picture ex- 
periences corrections when one moves through the phase 
diagram. Namely, in some areas of the phase diagram 
the splitting is extremely small, while in some others it 
may become so large that some multiplcts are pushed in 
the continuum and become unstable. The phase diagram 
presents different types of quantum critical points. At 
second order transition lines masses of certain particles 
vanish. Very close to the first order transition line addi- 
tional generations of kinks emerge. Strong interactions in 
some sectors may generate additional bound states (like 
breathers in the asymmetric charge sector). 

As we have mentioned many times throughout the text, 
one potential application of this theory is 'telephone num- 
ber' compound Sri4_ 2 .Ca x Cu2404i. This is undoubtedly 
a strongly correlated system. The measurements of low- 
frequency dielectric and optical response demonstrate ex- 
istence of a weakly pinned phason mode^^ which we 
identified as $c + ' mode. At the same time such probes as 
NMR 3 , inelastic neutron scattering 4 and ARPES 5 show 
gaps in all other parts of the spectrum. This is in agree- 
ment with existing theoretical understanding of the lad- 
der materials. 



The question of validity of the quasi-one-dimensional 
field theory description is decided by (i) comparison be- 
tween the values of the gaps and the bandwidth and by 
(ii) presence of essentially one-dimensional effects, such 
as different gap values for different channels. The field 
theory is valid when the gaps are small compared to the 
bandwidth. The gaps extracted from the ARPES mea- 
surements were obtained only for x = where the num- 
ber of holes is apparently rather small (though not zero, 
since according to Ref. |3 the gapless CDW mode exists 
at this concentration). The ARPES shows the single- 
electron gap ~ 0.3eV and the bandwidth ~ 1.2eV— . At 
the same time the neutron scattering (also available only 
for x = 0) gives the spin triplet gap ~ 32 meV with the 
bandwidth for the spin excitations ~ 200 meV. NMR 
measurements done in a broad range of Ca concentra- 
tions show that the spin gap decreases with doping and 
becomes less than 200K at x > 3. The absence of temper- 
ature saturation of the magnetic susceptibility indicates a 
crossover to the paramagnetic regime meaning that gaps 
for non-magnetic excitations also become smaller at these 
x. Therefore the gap/bandwidth ratios are sufficiently 
small for the field theory description to be valid. 

On the other hand, the optical conductivity measure- 
ments indicate a presence of essentially one-dimensional 
effects. They show a strong peak at radio frequen- 
cies (presumably coming from the phason mode) and 
a threshold at infrared frequencies (the so-called CDW 
gap) . If CDW state would form as a result of Fermi sur- 
face instability, as it happens in three-dimensional sys- 
tems with nested Fermi surfaces, this gap would be twice 
as large as the single-particle gap measured by ARPES. It 
would also coincide with the spin gap. However, accord- 
ing to Ref. 59, the values of the CDW gaps at x = 0, 3 and 
9 are 130 meV, 110 meV and 3 meV respectively which is 
either several times larger or much smaller than the spin 
gap. From this analysis we conclude that though de- 
tailed comparison between theory and experiment would 
be premature, the field theory description of telephone 
number compound is a reasonable approach ( as far as 
the system remains strongly one-dimensional which prob- 
ably corresponds to x < 9). 

We believe that the present experiments allow one to 
determine the part of the phase diagram where the tele- 
phone number compound is located. The recent X-ray 
measurements show that the holes may crystallize in a 
three-dimensional Wigner crystal3*Si (the use of the term 
is explained in^). Since the particle peaks sharpen be- 
low the transition, we take it as an indication that this 
is a Wigner crystal replacing the SCd phase, as de- 
scribed in Section IIIII Such crystal exists only if the 
electron-electron interaction has a long range tail. This 
points towards the standard Coulomb interaction as the 
primary agent of its formation. Such interaction may 
lead to relatively small values of 1 and consequently to 
new bound states. It is just possible, that one such bound 
state (a breather in the ' sector) appears as a sharp 
peak in the Raman scattering experiments^ 1 -. 
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APPENDIX A: BOSONIZATION 



We adopt the following notations: 



v,a y/2iva 



%/27ra 



(Al) 



where p = ±1 and ip and (p are bosonic fields with right 
and left chirality. The Klein factors satisfy anticommu- 
tation relations of the 0(4) Clifford algebra 

{Va,Vb} = Sab 

and can be chosen as follows^: 

v-iavi,=h n Pl v P i=K-i) (p+1)/2 - (A2) 

The chiral bosonic fields are decomposed into the normal 
modes as follows: 



■P<? 



(A3) 



with the same decomposition for tp. The bosonic field, 
<&, and its dual, O as usual are 

- + m - ©i ±) - - $3 ■ (A4) 

The Ising model order and disorder parameter fields 
a, fi are related to the bosonic fields as follows: 

cos[\/7r$a ^] = CTif72, sm[-\/7r3>y^] = fiifm, 
cos[V^Oi +) ] = ^i(T 2 , sin[V^Oi +) ] = 01//2, (A5) 

with the similar formulas relating $i ' , oi ' to 

0"3, M3) MO- 



APPENDIX B: SEMICLASSICAL ANALYSIS OF 
GENERALIZED 0(2n) GN MODELS 

In order to have a more intuitive picture of the effect 
of perturbations (|26|l and (|28|l on GN models we repeat 



the semiclassical analysis of Sec. for N and M even, 
N = 2n and M = 2m, following Ref. The advantage 
of this assumption lies in the fact the Majorana fermions 
can be bosonized in couples and then the potential of the 
0(2n) GN model Q takes the form 



V 



cos(\/47r0 i ) 



(Bl) 



It has two families of minima, 4n — y/nn and <f>i — 
\pH(n + 1/2) that are usually referred to as positive and 
negative vacua because they correspond to positive and 
negative values of ^ a V'kV'l = Si cos v^irfa- Fermionic 
excitations (when stable) correspond to configurations 
interpolating between minima that belong to the same 
family, for instance from a configuration (0,...,0) at 
x — ► — oo to (±y/n,0, . . . ,0), at x — ► +oo. On the 
other hand, kinks interpolate between minima of differ- 
ent families, for example, from (0, . . . , 0) for x — > — oo to 
(±y/n/2, . . . , ±y/iv/2) for x — > +oo. It is easy to see that, 
while there are N possible fermionic states, the number 
of kinks states is 2 n . Please note that since kinks interpo- 
late between positive and negative vacua the mean value 
of ipftipl changes sign along a kink configuration. 

From this picture it is clear that fermions can be con- 
sidered as bound states of kinks. If fact, for instance, 
an elementary excitation associated to the transition 
(0, . . . , 0) — > (y/iv, 0, . . . , 0) can be obtained also with 
two transitions (0, . . . , 0) -> (+^n/2, +y/n/2) -> 
(v^Tj 0, . . . , 0) , where each of the two jumps corresponds 
to a kink. The stability of fermions against the decay 
into a pair of kinks depends on N. 

Let us now first consider the effect of a simple pertur- 
bation of the form 



V 



^ COS V^TT(j)i 



(B2) 



Here the positive and negative minima are split and 
only the former remain absolute minima. This im- 
plies that kinks interpolating between different fami- 
lies are not present anymore and the model has only 
fermionic excitations. This is a simple example of 
confinement. If we now consider the model i|28[l . 
the effect of the perturbation on the excitations of 
the two GN models, originally decoupled, is somewhat 
similar. When the two models are decoupled each 
of them has the same families of minima described 
above. For instance, if we introduce <f>i and <Ji associ- 
ated to the x°s and £ & s respectively, configurations like 
(0i,...,^„;o-i,...,o-„) = (O,...,O;±0F/2,...,±V^/2) 
are absolute minima. Nevertheless this does not remain 
true in presence of a perturbation of the form 



7a. 6 < 



COS V 4:1V (j) a COS V 4lV(7 



(B3) 



(here we assume that also the perturbation is such that 
it can be written in a simple bosonic form). Again 
the kinks of the two decoupled GN models confine 
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and disappear from the spectrum. At the same time, 
together with fcrmions, new kinks with the quantum 
number of 0(N + M) GN appear. They interpo- 
lates between configurations like (0, . . . , 0; 0, . . . , 0) — > 
(±0F/2, . . . , ±Vtt/2; ±Vtt/2, . . . , ±y/n/2), that remain 
degenerate minima also of the perturbed theory. 
Fermionic excitations of the two models remain stable. 

Following the same arguments one can easily see that 
in presence of perturbations of the form l|26(l the O(N) 
kinks remain stable. In fact, the potential 

Sg a .b cos VATT<p a cos y/4Tr<pb (B4) 

does not lift the degeneracy of the minima of the unper- 
turbed GN model. 



particle and one unphysical at 27ri/3. Each of the S- 
matrices is crossing symmetric, including the scalar fac- 
tor£: £(0) = -£(m - 9). 

For the (3+1) model we have a similar structure: 

*«-(^«) «» 

where 



and S su ( 2 > is the S-matrix of SU(2) Thirring model soli- 
tons. 



APPENDIX C: THE EXACT SOLUTION OF 
(5+1) AND (3+1) MODELS 

In this appendix we describe the S-matrix of two mod- 
els, that we call (5+1) and (3+l)-model, with 0(5) x Z2 
and 0(3) x Z2 symmetry respectively. The l atter was in- 
troduced and studied in detail in Ref.s l42l43L The central 
charge in the UV computed with the Thermodynamics 
Bethe Ansatz is C5+1 = 3 and c^ + i = 2. The relationship 
with the models of interest in the paper is discussed in 
Sec. El 

The spectrum of the (5+l)-model consists of a singlet 
Majorana fermion with mass mo, a quintet of Majorana 
fermions with masses y/3M and a quartet of kinks = 
antikinks (spinor particles) with mass M. The S-matrix 
for the (5+1) model is 



5., 



S VS S SS £(5^ 

-I -1 



(CI) 



where 



a 30 



and S vv , S vs and S ss are the 0(5) S-matrices of vector 
and spinor particles. The above S matrices were found 
in Ref. |38| (for the vector particles) and |39| (for kinks) . 
The spinor S-matrix has the following form: 



S ss (9) = f(8) 



p 

1 asyn 



iyr/3 



i?r/3 



P, 



+ m 



Pn 



,(C2) 



where Po,P v ,P a sym represent projectors onto singlet, 
vector and antisymmetric tensor representations and 



/(<?) = 



rd + f r |-f 



H r i + f 



r(i-i)r(i + fe)r(f + t)r(|-i) 



The S-matrix has one physical pole on the physical strip 
(0 < ^sm8 < 7i") at 8 — in/3 corresponding to the vector 



APPENDIX D: SOLUTION OF THE 
LINEARIZED RG EQUATIONS 

The RG equations JSJ can be simplified for g a = g+Sg a 
with g/Sg a <C 1. Defining for simplicity of notations 
Sg a = x a we find that, to first order in x Eq.s JSJ take 
the form 



= -Dx 



where 



D = 



2.9 



/ 1 3 0\ 

' 1/2 1/2 3/2 3/2 ' 

1111 

1/2 1/2 1/2 3/2 1 

V 1 21/ 



(Dl) 



(D2) 



and g is solution of g = — 4g 2 and has the form 



9(t) = 



1 



At 



go 1 



(D3) 



It is convenient to introduce y 

x = Ty, (D4) 
that satisfies the equation 

y = -T- 1 DTy2g(t), y(0) = T -1 x(0) = A a , (D5) 
where T 



T = 



/ -3 -6 -1 2 1 \ 
'3 3-1-11' 
-2 -2 0-2 1 
10 11 
V 1 1 1/ 



(D6) 



is chosen such that T~ 1 DT = diag(Ai, A2, A3, A4, A5), 
with As being the eigenvalues of (ID2|) . Ai = A2 = — 1, 
A3 = A4 = 1, A5 = 4. Since the eigenvectors of degener- 
ate eigenvalues are linearly independent all solutions of 
(ID5II have the form 



Va = A a exp[-\J(t)] 



(D7) 
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where where the relationship between t and e was used. The 

j solution of (|D1|) can be obtained inverting (|L)4|i . 



From this it follows that 



1 



Aa/2 



!la = A a [^-) (\og(e/m) + iy xj2 , (D9i 
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